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A general geometrical structure of the entanglement entropy for spatial partition of a relativistic 
QFT system is established by using methods of the effective gravity action and the spectral geometry. 
A special attention is payed to the subleading terms in the entropy in different dimensions and to 
behaviour in different states. It is conjectured, on the base of relation between the entropy and the 
action, that in a fundamental theory the ground state entanglement entropy per unit area equals 
l/(4Gjv), where Gjv is the Newton constant in the low-energy gravity sector of the theory. The 
conjecture opens a new avenue in analogue gravity models. For instance, in higher dimensional 
condensed matter systems, which near a critical point are described by relativistic QFT's, the 
entanglement entropy density defines an effective gravitational coupling. By studying the properties 
of this constant one can get new insights in quantum gravity phenomena, such as the universality 
of the low-energy physics, the renormalization group behavior of Giv, the statistical meaning of the 
Bekenstein-Hawking entropy. 



PACS numbers: 04.60.-m, 03.70.-fk, 03.65.Ud, 05.50.-|-q 



I. FUNDAMENTAL ENTANGLEMENT 
A. General properties 

Quantum entanglement is an important physical phe- 
nomenon in which the quantum states of several objects 
cannot be described independently, even if the objects 
are spatially separated. Quantum entanglement is used 
in different research areas P]. In quantum information 
theory entangled states are a valuable source of process- 
ing the information. Quantum entanglement also plays 
an important role in properties of strongly correlated 
many-body systems and in collective phenomena such as 
quantum phase transitions. In quantum gravity the en- 
tanglement may be a key for understanding the mystery 
of the black hole entropy. 

The entanglement can be quantified by an entropy. 
One can define it as the measure of the information about 
quantum states which is lost when these states cannot be 
observed. In many-body systems, which are the subject 
of the present work, "observable" and " unobservable" 
states can be located in different regions. Consider, for 
instance, a lattice of spins being in a quantum state char- 
acterized by a density matrix p. Suppose that the lattice 
is divided into regions A and B with a common bound- 
ary B. The entanglement between the two regions can 
be described by the reduced density matrix 

pB = Tr^/5 , 

where the trace is taken over the states of spin operators 
at the lattice sites in the region A. The entanglement 



entropy in the region B is defined as the von Neumann 
entropy 



Sb = -TtbPb Inps 



(1) 



Analogously one can define the entanglement entropy 5*^ 
in the region A by tracing the density matrix over the 
states located in the region B. If the system is in a pure 
state 1-0), i.e. p_= it is not difficult to show that 

Sa = Sb: see [2. 

It has been observed long ago 0]j[3 (see also recent 
works 0, 0), that in lattice models representing a dis- 
crete version of a relativistic quantum field theory (QFT) 
the ground state entanglement entropy in the leading or- 
der is proportional to the area A of B, 



S — Sf 
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(2) 



Here g is the lattice spacing and D is the number of 
spacetimes dimensions (it is assumed that D > 2). This 
geometrical property of S follows from the fact that field 
excitations in different regions are correlated across the 
boundary. 



B. Black hole entropy 

If g is identified with the Planck length Ipi, (O looks 
very similar to the Bekenstein-Hawking entropy of a 
black hole 
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Here Ah is the area surface of the black hole horizon and 
Gn is the Newton constant, G^v = ^p;"^- The similar- 
ity of the two entropies suggests that S^^ may be an 
entanglement entropy, as was first pointed out in |^,[^. 
An external observer, who is at rest with respect to a 
black hole, never sees what happens inside the horizon. 
Because the observer perceives the vacuum excitations in 
a mixed state it is natural to assume that S^^ measures 
the loss of the information inside the horizon. 

The reduced density matrix for a black hole is thermal 
. Thus, the entanglement entropy coincides with the 
entropy of a thermal atmosphere around the black hole 
horizon (see for a review of different interpretations 
of S^^ in terms of the entropy of quantum excitations) . 
However, relation between S^^ and the entropy of en- 
tanglement (or thermal entropy) is non-trivial because S 
is divergent. In the lattice regularization ^ diverges in 
the continuum limit when g ^ 0. It was pointed out in 
0, that the divergences in S are the standard ultra- 
violet divergences. If the entropy of thermal atmosphere 
is considered as a quantum correction to S^^ the diver- 
gences can be eliminated in the course of renormalization 
of the gravitational coupling. 

A deep connection between the divergences in S and 
the divergences in the quantum effective action has been 
established latter for different field species at the one- 
loop level 7]. This connection is very important. Sup- 
pose, following the original idea by Sakharov that 
the Einstein action is entirely induced by quantum ef- 
fects of some underlying fundamental theory. Because 
the Newton constant in this case has a pure quantum 
origin, so does the entropy S^^ . The hypothesis is that 
the entropy of a black hole in such a theory is an en- 
tanglement entropy of fundamental quantum excitations 
which induce low-energy gravity dynamics [TH-|l3j. 

The mechanism of generation of the black hole en- 
tropy in induced gravity may be quite general and, 
in particular, it can be realized in string theory |l4j . 
In string theory the low-energy effective gravity action 
appears from the tree- level diagrams of closed strings. 
From the point of view of open strings these are one- 
loop diagrams. Thus, the gravity action as well as the 
Bekenstein-Hawking entropy can be interpreted as a pure 
loop effect. Another argument in favour of entanglement 
interpretation of black hole entropy comes from analysis 
of string theory on asymptotically anti-de Sitter back- 
grounds and its duality to a conformal field theory (CFT) 

MM 



C. Entanglement entropy in fundamental theories 

The purpose of the present work is to look at these 
results from a different perspective. Suppose the under- 
lying fundamental theory of quantum gravity is known. 
This can be a string-I?-brane theory or some other theory 
which correctly describes the observable particle physics 
and gravitational effects. Consider a low-energy limit of 



this theory and the entanglement entropy S in the ground 
state for partition of the system in Minkowsky spacetime 
by a plane. It is assumed that S measures entanglement 
of all genuine microscopical degrees of freedom of the 
fundamental theory (which may not coincide with the 
low-energy fields). We call such quantity the fundamen- 
tal entanglement entropy. 

Because the space is infinite it is convenient to discuss 
the density of S per unit area. Let us denote this quantity 
by S. We make a conjecture that the density of the 
fundamental entanglement entropy is finite and is given 
precisely by the relation 

'^ik ■ 

where Gn is the Newton constant in the low-energy grav- 
ity sector of the theory. The conjecture is based on the 
above mentioned arguments: 

i) in a fundamental theory the cutoff parameter g in 
(0) should be of the order of the Planck length; 

ii) there are evidences that the Bekenstein-Hawking 
entropy ^ is the entanglement entropy for partition of 
the system by the horizon. 

In what follows we discuss a more formal argument in 
favour of (@J). It is based on the observation that the 
black hole entropy and the entropy of entanglement are 
derived in a unique way from the effective gravity action. 



Relation Q has a number of interesting consequences. 
One of them is that the definition of Gat becomes not only 
the subject of gravitational physics. The flat space entan- 
glement is an alternative source of information about the 
Newton constant. (In a certain sense the situation with 
different definitions of Gat can be compared with "grav- 
itational" and the "inertional" notions of the mass.) 

In this paper we pay attention to the fact that re- 
lation I0J) can be used to study the properties of the 
gravitational coupling by carrying the computations of 
the entanglement entropy in condensed matter systems. 
There are higher-dimensional condensed matter systems 
which near a critical point corresponding to a second or- 
der phase transition are described by relativistic QFT's 
which contain massive fields. The class of Ising models 
is one of such examples. The QFT's in these systems are 
analogous to low-energy limit of the fundamental the- 
ory. The relativistic covariance guarantees that the en- 
tanglement entropy density in the near-critical regime de- 
termines an effective gravitational coupling by the same 
equation as Eq. Q). The lattice spacing in the under- 
lying theory serves as a cutoff. Hence the entanglement 
entropy is finite and can be calculated by using analytical 
and numerical methods. This enables one to pose differ- 
ent questions about properties of Gjv and try to find the 
answers in gravity analogs. 
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D. Outline of the paper 

The rest part of the paper is organized as follows. Sec- 
tion II is devoted to the entanglement entropy in rel- 
ativistic QFT's. The Section starts with discussion of 
one-dimensional spin models which serve as an illustra- 
tion. We then establish connection between the entan- 
glement entropy S in relativistic QFT's and the effective 
gravity action, and use methods of the spectral geometry 
for derivation of S. By studying the theories in a cubic 
region we find a geometric form of S in different dimen- 
sions, at zero and high temperatures. We pay a special 
attention at boundary effects. Our result for the ground 
state entanglement in four-dimensional space-times looks 
as follows: 
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where A is the area of the separating surface B and C 
is the length of the boundary of B. The last term in S 
has a topological form. It depends on the number C of 
sharp corners (vertexes) of the boundary of B and does 
not change under smooth deformations oiB. Calculation 
of the subleading terms in is the new result. 

Relation of S to the Bekenstein-Hawking entropy in 
induced gravity models in four dimensions is discussed in 
Section III. We also discuss the concept of gravity analogs 
and possible models which can be used to study proper- 
ties of the gravitational coupling on the base of relation 
Q . There are several topics which can be pursued by us- 
ing simulations in the new class of gravity analogs. Some 
of them such as the renormalization group behavior of 
Gjv, the universality of the low-energy physics and the 
statistical meaning of the Bekenstein-Hawking entropy 
are listed in Section IV. 



II. ENTANGLEMENT IN RELATIVISTIC 
FIELD THEORIES 

A. One-dimensional spin chains 

We begin with discussion of the properties of the en- 
tanglement entropy in one-dimensional spin chains. Con- 
sider as a simplest example the entanglement entropy for 
a block of spins in the Ising model. The Hamiltonian of 
the model is 



N 



(6) 



where N is the number of spins and af , are the Pauli 
matrices. Parameter A is the strength of external mag- 
netic field. At zero temperature the Ising chain has a 
second order phase transition at the critical value A = 1. 

As is known, there is a unitary transformation (the 
Jordan- Wigner transformation) which maps the spin 



models like (jHl to one-dimensional models with two spin- 
less fermions. This transformation can be used to diago- 
nalize the Hamiltonian. 

The entanglement entropy can be investigated for the 
ground state when the chain is separated into two blocks 
of contiguous spins of equal sizes N/2. Suppose that A is 
fixed and N varies {N is much larger than the correlation 
length). The behaviour of S has been studied in two 
regimes In the off-critical regime, A ^ 1 and |A — 

1| ^ 1, the entropy at large N reaches the saturation 
value 



S{N, A) 



log2|A-l| 



(7) 



In the critical regime, A = 1, the entropy does not reach 
the saturation and behaves at large N as 



SiN, A)~ ilog2Af/2 
6 



(8) 



The explanation of these results is based on the fact that 
in the continuous limit, iV — > oo, the Ising model near 
the critical point corresponds to a quantum field theory 
with two fcrmion fields whose mass m (or the inverse 
correlation length ^ ) is monotonically related to | A — 1 1 . 
In the critical regime m vanishes and one has a conformal 
field theory with two massless fermions each having the 
central charge 1/2. 

Formulae ((HJ are quite general. They have been 
verified in other exactly solvable one-dimensional statisti- 
cal models (XXZ, XY and other models) which generalize 
the Izing chain see and a review of results in 

[23|. To discuss the continuum limit in lattice models it 
is convenient to introduce the lattice spacing g. Then the 
number of spins N defines the size of the system L = Ng. 
In the systems separated into two equal parts the entropy 
at the critical points behaves as 



c L 
-5" - T In — 
6 2g 



(9) 



where c is a total central charge of the corresponding ef- 
fective two-dimensional conformal theory. Near the crit- 
ical point the entropy behaves as 



6 gm 



(10) 



where m is the corresponding inverse correlation length. 
From the point of view of QFT g is the ultraviolet cutoff 
parameter. 



B. Entanglement entropy and effective gravity 
action 

There are different ways to derive ©, (fTUIl by QFT 
methods. Our purpose is to show how it can be done by 
using the effective action approach. We first start with 
a theory in a flat spacetimc with arbitrary number of 
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dimensions D and then focus on theories in D = 2,3, 
and 4. 

The entanglement entropy for a lattice system dis- 
cussed in Section LA can be rewritten as 



d 

Sb ^ - hm tt-Ttb Pb 
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(11) 



Suppose that near a critical point the system is equiva- 
lent to a QFT with field variables </> whose dynamics is 
determined by the action /_b [(/)]. The density matrix pB 
in the configuration representation depends on variables 
(j)B in the region B. We consider the system at a finite 
temperature T. The result for the ground state entan- 
glement can be obtained in the limit T — *■ 0. The matrix 
elements of can be described in terms of the Euclidean 
path integral 



Pb\<Pb) =AA- 



[D4>] 



(12) 



where A/" is a normalization coefficient introduced to sat- 
isfy the condition Tr^ pB = 1- The classical action I 
is defined on a Euclidean space with Euclidean time t 
compactified on a circle of length T^^. To be more 
specific suppose that the partition of the space is done 
orthogonally to one of the spatial coordinates, say x^. 
Then regions A and B can be defined, respectively, as 
-~L/2 < < and < < L/2, the separating sur- 
face is located at x^ = 0. Relative sizes of A and B can 
be arbitrary in general. The integration in l|12|l goes over 
field configurations (j) with the boundary conditions 



.,x 



0) = 0B(a;\...,a;^-^) 



,x^-\t = T-') = M^\ 



where a;^ takes values in the interval (0, L/2). 

Let us denote by the space where the field config- 
urations (j) are given. It is easy to see how AA'-^ looks in 
two dimensions. It is a cylinder (see the upper picture 
on Fig. ^ with coordinates {x^,t) and a cut parallel 
to its axis, the length of the cut is equal to the length 
of the interval B. The field takes values (t>B, (t>'B the 
lower and upper sides of the cut, respectively. In higher 
dimensions the cut is a (Z? — 2)-dimensional hyperplane. 

If the parameter n in Hll|) is positive and integer the 
matrix elements of the operator p^ can be represented 
by integral \i2l where field variables are defined on a 
space M.'n which is obtained by gluing n copies of M.'i 
along the cuts. M'^ still has a single cut which disappears 
when one takes the trace, Trs pg . We denote the space 
obtained as a result of this procedure as A^„. The lower 
picture on Fig. ^ shows in two dimensions. 

Ain is locally flat but has a non-trivial topology. Its 
boundary in two dimensions consists of n circles of the 
length T^^ each and a circle of the length nT~^ . The 
important property is that Ain has a singularity at a 
hyper-surface where all n cuts meet. For n>2 each point 



1/T 




3/T 



FIG. 1: The upper picture shows M'^ in two dimensions. It is 
a cylinder with the circumference length and a cut along 
the axis. The space is schematically drawn on the lower 
picture. It is obtained by gluing along the cuts of 3 copies of 
M'l. The circumference length of the right boundary of A^s 
is 3r~^. The cuts meet at the point X which is a conical 
singularity. 



on this hyper-surface is the conical singularity because a 
unit circle around it has the circumference length equal 
27rrt. 

The singular behaviour of Ain is easy to see in two 
dimensions. A4„ has the same topology as a disk with n 
holes. Thus, its Euler characteristic. 
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47r 
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equals 1 — n. The boundary 9A^„ of A^„ consists of n + 1 
circles, each of which has a vanishing extrinsic curvature 
A; = 0. Therefore, the integral of the scalar curvature R 
for M-n should be non-trivial for n > 1, 



/ R = 47r(l - n) 



IM 

The non-zero value of the integral is ensured by the con- 
ical singularity. If this point is located at x = the 
curvature is the distribution R = 47r(l — n)(5(^^(a;), see, 
e.g. M- , 

Let us introduce the parameter (3 — 27m and rewrite 
in another form following from (|12|l . 



(13) 



r(/3) is the effective action defined by the path integral 



-r(/3) _ 



(14) 



where field variables are given on A4„. i^; [</),/?] is the 
classical action on A4„. The normalization of the integral 
([Tl|l is fixed by the condition r(27r) = - In TV. 
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The operation with the parameter /3 in H13() should 
be understood in the following way: one first computes 
r(/3) for 13 = 27rn, and then replaces /? with a continuous 
parameter. This can be done even if A^„ itself cannot be 
defined at arbitrary /3. 

If one defines a "free energy" F(/3) = r(/3)//3 and in- 
terpretes (3 as an inverse temperature, equation for- 
mally coincides with the definition of the entropy in sta- 
tistical mechanics. Note, that (3~^ is a "geometrical tem- 
perature" which should not be confused with the physical 
temperature T. 



C. Spectral geometry on Mn 

The expression of the entanglement entropy in terms 
of the effective action has a number of advantages. It 
enables one to reformulate the problem on the geomet- 
rical language and use powerful methods of the spectral 
geometry to study the form of the entropy. 

Consider a Laplace operator A = — on 7W„ and the 
trace of its heat kernel, Tr e~*^, where Hs a positive pa- 
rameter. The spectral geometry relates the asymptotic of 
the trace at small t to integral geometrical characteristics 

of Mn, 



Tr e 



1 



(47rt)^/2 



(15) 



where n takes non-negative integer and half-integer val- 
ues. The coefficients a„ are expressed in terms of the 
powers of the Riemann tensor and its derivatives, and 
include boundary terms (see, e.g., H^). The coefficients 
also depend on the singularities of the background man- 
ifold such as jumps of the curvature, sharp corners, and 
so on. The effect of conical singularities appears starting 
with the coefficient Oi . For scalar and spinor Laplacians 
in two dimensions 



ai 





(2ti 










27r J 
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(16) 



For spaces A^„ like those shown on Fig. ^the extrinsic 
curvature for each element of the boundary is zero and oi 
is determined only by the conical singularity. Other non- 
vanishing coefficients in the asymptotic expansion for the 
heat kernel on A^„ in D = 2 are 



ao = lvol{Mn) = ljr = 



(17) 



ai/2 



OMr, 



(18) 

For scalars and spinors I = 1 and 1 = 2, respectively. The 
numerical coefficients b'^^^ ^^'^ ^i/r> ~ ~^t/2 correspond, 
respectively, to the Dirichlet and Neumann conditions 
on the boundaries of Al„. For the scalar Laplacian = 



— 1/4 and b^ = —1/2 for the spinor Laplacian. The non- 
trivial dependence on f3 in the first term in the r.h.s. of 
()16|l is the direct consequence of the conical singularity. 
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FIG. 2: This figure shows division of a cube by a plane B. 
The trace is taken over the states of a QFT in one of the halts 
of the cube. 



Consider now QFT in a (_D— l)-dimensional cube with 
the Dirichlet or Neumann boundary conditions on its 
faces. Let L be the length of the edge of the cube. Sup- 
pose we are interested in the entanglement entropy asso- 
ciated with the partition of the cube by a plane B which 
is orthogonal to one of the faces and divides the cube 
into two equal parts, as is shown on Fig. [21 Let us de- 
note by A^i^^ a Z)-dimensional manifold which appears 
in calculation of the entanglement entropy in terms of 
the effective action along the lines explained in Section 
ILB. 

In two dimensions Ain'' is shown on Fig. ^ If _D > 2 
the manifold M^^^ is the space product of Mn ^ and 
D — 2 intervals of the length L. For a scalar Laplace 
operator A^"^^ on Aii!^^ one can write 



Tr 



D-2 



Tr e 



(19) 



where KL{t) is the trace of the heat kernel of the Laplace 
operator on the interval HJ. The following asymptotics 
hold at small t up to exponentially small terms 



Tr e 



-tA(=' 



4^("° 



ai/2V^ + fliO 



1 



where a„ are given in p6|l - (|18|l . Therefore, 

,-tA(°> 1 



Tr e 



^ ' n=0 



D/2 



(20) 



(21) 



(22) 



where coefficients al^'' can be found from (|20|1 . l(2U- One 
can check with the help of (|16fl - (|18|l that the first coeffi- 
cients take the form 



2-kT' 



L 



D-l 



= vol(A^(^)) , 



(23) 
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(24) 
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Li-4 



(26) 




FIG. 3: The manifold A^g is shown for computing the entan- 
glement entropy Sb at a finite temperature when the region 
A consists of two internal disjoint intervals. M'^P is obtained 
by gluing 4 copies of the manifold shown on Fig. -M^^^ has 
4 conical singularities at branch points. 



(£)) _ TT / 27r (3 



4^&±,)4^2= 



ttT" 



(I?-l)! 



Li-5 



(27) 



4!(i:'- 5)!^ 

The heat kernel coefficients for the spinor Laplacian have 
a similar structure and can be found in analogous way. 

Although our computations have been done in a sim- 
ple setting, the results can be easily generalized when 
M']^^ is an arbitrary curved manifold. We are particu- 
larly interested in the geometrical form of contributions 
to the heat kernel coefficients from the conical singular- 
ities because they yield a non-trivial value of the entan- 
glement entropy. These contributions depend on the fac- 
tor ((27r//3) — {P/2tt)) and appear in the coefficients a„ 
with n > 1. The geometrical structures related to con- 
ical singularities in the coefficients Oi, 03/2, and 02 are, 
respectively, 



L 



D-2 



2{D-2)L 
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(28) 



(29) 



D. Effective action and the entropy 

Let us show now how the geometrical structure of the 
entanglement entropy can be established by using the 
effective action. We consider, as an example, a free scalar 
field theory in a cube divided into two parts, Fig. |21 The 
effective action for a massive scalar field can be defined 
as follows 



1 



-tA-tm^ 



r = -lndet((A + m^)^-^) = — / — Tr e 
2 2 Jg2 t 

(31) 

where /i is a dimensional parameter. The right hand 
side of gives the effective action in the Schwinger- 
De Witt representation. The parameter g is & ultraviolet 
cutoff which has a dimensionality of a length. 

It is assumed that the correlation length m"^ of the 
theory is much larger than g. Let us also assume that 
TO~^ is much smaller than the size of the system L. In 
this case one can use for F an approximation which can be 
obtained when the heat kernel is replaced by its asymp- 
totic if^ . 

Here and in what follows we use instead Sb the nota- 
tion S because in the considered examples the regions A 
and B are identical. 



2{D -2){D -3)L 



D-i 



= C 



dB 



(30) 



The separating surface S is a (D — 2)-dimensional man- 
ifold. In general, it may have a boundary dB and the 
boundary may consist of different hypersurfaces where 
the curvature of dB has jumps. Integrals in the right 
hand sides of (|29l) and (|3U|I are taken, respectively, over 
the boundary dB and codimension 1 singular hypersur- 
faces inside dB. In the considered example of the theory 
in the hypercube, K is a square if _D = 4 (see Fig. EJ. It 
has the area A = \i the edge length is L. Its boundary 
dB has the length 4L and consists of 4 singular points, 
the vertexes of the square, li B \s a. D — 2-dimensional 
hypercube, 2{D ~ 2) is the number of its "faces" and 
2(D - 2)(D - 3) is the number of its "edges". 



1. D=2 
In two-dimensional theories 
F^-A 



(2) -2 

% 'q 



2af/2i?" 



, , 2 (2) (2)n , / 2 2\ / (2) 2 (2) 1—^ 

+ (m Uq' — a\ ) m(m g ) — (Oq 'm — a^jr^m^/'n) 

;(32) 

In this expression we drop all terms which vanish in the 
limit — > 0. To find the entropy one can apply H13|l to 
(|^ and use It yields 



(33) 
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If Q is identified with the lattice spacing, this resuh coin- 
cides with computations in spin chains in the near critical 
regime, see (|1(JII . The central charge in the given case is 
c — \. Formula (|33|l is obtained for the partition of the 
interval into two parts. It can be easily generalized for 
more complicated cases when the regions A and B con- 
sist of several disjoint intervals. The effective action in 

(2) 

this case is defined on a 2D manifold Mn with some 
number of handles. Figure O shows M'"^'^ for tracing the 
degrees of freedom in two disjoint internal intervals. As 
is easy to understand the number of conical singularities 
Mc of equals the number of internal points of the 

intervals. Each conical singularity yields an independent 
contribution to the entropy. Thus, in general case the 
result in the right hand side of H33|l has to be multiplied 
hyNc |2|. 



2. D=3 



In three-dimensional theories 



r — 



167r3/2 



(3) -2 



2{a\ 



(3) 



)g ^ — (13/2 ~ ln(?7i^p^)-f 



^ (3) 3 (3) 2 r- (3) 
——-=aQ 'm — aiLm — y'na\ m 



and the entanglement entropy takes the form 



(34) 



1 



1 



(35) 

where we used (|28|l . (12911 . For the theory in the square 
A = L, £ = 2. The term in (|35|) which is proportional 
to C has the following property: it does not change un- 
der continuous deformation of the separation surface B. 
This property is in agreement with a more general be- 
haviour of the entanglement entropy in three-dimensional 
theories which has a topological part arising from short 
wavelength modes localized near the boundary of B, see 



■^ 4 4 , 

4 ^ 



(4) 3 (4) 2 



,(4) 



(36) 

Equations (|23(l - (|30|l give the following expression for the 
entanglement entropy 



s. 



487r 



+ m^{ln{m^g^) - 1) ] A- 



/2 



'J 'Km 



C — — ln(m^p^ 



(37) 



For a theory in the cube, C = 4L, C = 4. In general, C 
is the number of singular points of the boundary dB of 
the separation surface. The term in H37I) which is propor- 
tional to C does not change under the continuous defor- 
mation of dB. Therefore, it is a topological term similar 
to the subleading term in the entropy in three dimen- 



Expressions (|35|) . (|37|l hold for relativistic theories. In 
general, dispersion relations in condensed matter systems 
are not Lorentz invariant and are not linear. As a result, 
the leading term in the ground state entanglement en- 
tropy may have a more complicated dependence on the 
area of the partition surface A than in (|35|l . H37() . In par- 
ticular, an additional logarithmic coefficient pointed out 
in H^-ll^ may appear in the scaling of the entropy. 



E. Ground state entanglement at critical points 

So far we assumed that the system has a non- vanishing 
correlation length which is small compared to the size 
of the system. At critical points the correlation length 
becomes infinite and approximations considered above do 
not hold. In general, the effective action in this regime is 
non-local and its computation is non-trivial. 

The effective action F is the sum 

r = rH,. + F_ . 



3. D=4 
In four-dimensional theories 



1 



327r2 



1 



(4). 



2«o 9 



(4) 
h/2f 



(4^) - m^aj^' + —4")) In(m2g2) 



Fdiv diverges in the limit of vanishing cutoff g while the 
" renormalized" part F^^n is finite at gi = 0. The structure 
of divergences F^iv can be found from (|22Jl, (|S3|l , (EHl • To 
obtain F^iv in the massless theory one has to leave in (|32|l , 
only terms depending on g and take the hmit 

m = 0. 

For critical systems whose size is parametrized by a 
single parameter the renormalized part can be found by 
using anomalous scaling properties of F^on. Let us con- 
sider as before a free scalar field theory in a cube of the 
edge size L. The effective action is F = i Indet A/i^^. If 
one changes the operator A to A' = e^'^A, where cr is a 
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constant, the renormalized part of the action for A' can 
be written as 



T' 



r,., 



=_D/2 



(47r)^/2 



(38) 



This result can be easily obtained in the dimensional reg- 
ularization The last term in the right hand side of 
this equation is known as conformal anomaly. In the con- 
sidered case the action is a function r„„(/3, L, T, /x). By 
using lp!H)l one can write 

r_(/3,L,r,/i) = 



TUP,Lo,{LT)/Lo,^l) 



In- 



L 



(39) 



(47r)^/2 Lo ■ 
Analogously, the entanglement entropy can be written as 



1 

48^ 



\n{L/g) C 



(46) 

where we omitted constants which do not depend on the 
parameters of the system. By comparing (I44(l - (|46|l and 
(123 1 (123 7 one can conclude that the size of the 

system L in the massless theory plays the role of the 
infrared cutoff analogous to the inverse mass. The result 
for the entropy in two dimensions is in agreement with 
the entanglement entropy in the critical spin chains, see 
@ and the earlier computations of the entropy in .30] . 

Although formulae H44|l -H46 () are obtained for the the- 
ory of a free field, they should preserve their form in the 
presence of interactions. The logarithmic terms related 
to the conformal anomaly should not receive corrections. 



F. High temperature limit 



(40) 



where S^i, and S',„„ are divergent and finite parts which 
are calculated, respectively, from r^iv and r^^,,, by using 
equation (|13(l . According to (|39|l . the renormalized part 
can be written as 



SUL,T,^i) = 5„„(Lo,(ir)/Lo,M) + 



(D) 
^n=D/2 



In- 



(47r)^/2 Lo 



"-"-;is. (4-')°" 



^n=D/2 



(D) 



D/2 



(41) 



(42) 



The part of the entropy S'„„(Lo, (Lr)/Loi m) is deter- 
mined by r„„(/3, Lo, iLT)/Lo,fi). 

To get the entropy in the ground state one has to 
go to the limit T ^ 0. It is known from the lat- 
tice computations that the entanglement entropy is fi- 
nite at zero temperature. This property implies that 
Srei, (Lo-, (LT) / Lq, fx) in the zero temperature limit should 
be a finite constant, C{Lq, fi), which does not depend on 
the size L. Therefore, 



5_(L,T = 0,m) 



^ln^ + C(Lo,.) . (43) 



By using equations (1221, (123, 1211), (EH and 

(|43|l one finds the ground state entanglement entropy of a 
massless scalar field for a theory in a cube in dimensions 
D — 2, 3, 4, respectively, 
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HL/g) 



(44) 



The asymptotic of the entanglement entropy can be 
also found in the opposite limit, TL ^ 1, T ^ m. In 
this limit the main contribution comes from the modes 
with the wave length of the order of T^^. The theory 
in this case is extensive. Consider first the theory on 
an interval L divided into two parts. The corresponding 
background manifold Ain ^ is shown on Fig. ^ As before, 

(2) 

let rro„(/3, L, T, n) be a renormalized action on Adii {(3 = 
27rn). Let r(/3o, L/2, T, /i), f3o — 2tt, be the action on 

(2) 

the cylinder of the length L/2. The space Mii can be 
obtained by gluing n cylinders of length L/2 and the 
circumference length with the cylinder of length L/2 
and the circumference length nT~^. Because the theory 
in the high temperature limit is extensive the action can 
be written as 

r_(/3,L,T,Ai)~ 



nr(/3o, L/2, T, fi) + r(/3o, L/2, T/n, /i) 



(47) 



This relation also holds in higher dimensions. It is not 
difficult to show by knowing the free energy of a one- 
dimensional gas on an interval that at high temperatures 

f(/3o,L/2,T,Ai)^-5(L/2)r . 

6 

(The effective action is a free energy divided by temper- 
ature.) Therefore, 



r_(/3,L,r,//) 



12V^ 



26^^/2 ln(L/g) C 



(45) 



At /3 = /3o this result coincides with the effective action 
on the cylinder of the length L. The divergent part of 
the effective action in this limit can be neglected and the 
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entropy can be determined from (|48|l with the help of 
(113, 



S{L,T,y)^ -{L/2)T 



(49) 



The entanglement entropy in the high temperature limit 
coincides with the usual statistical entropy of the system 
of the length L/2. A similar result was found in j24j . 

In three and four-dimensional spacetimes the effec- 
tive action for a scalar theory in a volume V at high 
temperatures behaves, respectively, as — C(3)yr^/(27r), 
-7r2T/TV(90), see From one then finds high 

temperature entanglement entropy for partition of a 
square or a cube into two parts 



3C(3) 
2tt 



VT^ 



SiL,T,fi)^^VT^ 



(50) 



(51) 



Here ({x) is the Riemann zeta function and V = 
In general case, V in (|50|l . H51() is the volume of that re- 
gion where the entanglement entropy is calculated. Note 
that in this regime the entanglement entropy coincides 
with the usual statistical entropy in the volume V, in 
agreement with extensive properties of the system at high 
temperatures. 



III. FUNDAMENTAL ENTANGLEMENT AND 
INDUCED GRAVITY 

A. Induced Newton coupling 

According to Sakharov's idea the classical 

Einstein-Hilbert gravitational action can be entirely in- 
duced by quantum effects of matter fields on a curved 
background. As an example, consider a model of a free 
scalar field of the mass m with an ultraviolet cutoff g. 
The model is defined on a curved manifold A4 with met- 
ric g^i,. To keep connection with finite-temperature theo- 
ries we assume that 5^,^ has the Euclidean signature. For 
simplicity we also assume that is a compact closed 
manifold. The effective action of the field can be found 
by using (|31|l . The heat kernel operator now has the 
following asymptotic (see e.g. [2^ ) 



Tr e 



-tA 



1 



2 Z-^ ' 



JM D JM 



(52) 



(53) 



where g = det^^jy, R is the scalar curvature oi M., n 
takes non-negative integer values. The higher coefficients 



are integrals of powers of the Riemann tensor and its 
derivatives. 

We consider theory in 13 = 4. If the curvature radius of 
Ad is much larger than m^^ on can proceed as in Section 
II. D and get 



1 



327r2 



(54) 



In the last line all terms related to Ag^-* were omitted. 
The right hand side of 1)54(1 coincides with the classical 
Einstein-Hilbert action where G and A are the induced 
Newton and cosmological constants, correspondingly. 



A 
G 



1 

G ' 
1 

"8^ 



1 

12^ 



1 



+ m^{\n{m^Q^) - 1) 



2m^ 



+ 771^ ( \ii{n?Q^) - 



(55) 



(56) 



The constant G has the correct order of magnitude if the 
cutoff parameter g is of the order of the Planck length. 

By comparing (|^ with (p!7)l we conclude that, in ac- 
cord with the hypothesis |0J, the induced Newton con- 
stant in this simple model and the entanglement entropy 
of the system in the flat space-time in a cube are related 
as 



1 

G 



4 lim — 

L^oo A 



(57) 



Here A = L?' is the area of the separation surface. The 
result holds when the size of the cube is sufficiently large. 

The reason why takes place is because the conical 
space possesses a curvature concentrated at the tip. Let 
us consider the heat kernel coefficient ai^\ Eq. 
the theory on A^l^'' and compare it with 
If (3 in H25|) is close to 27r then 



51, for 
in inSJl. 



AD) 



i2(2.- 



-P)A + h.i. 



i?c-|-b.t. 



(58) 



where Rc = 2(27r - I3)5'^{x) is the "curvature" of M\ 
concentrated at conical singularities on the hypersurface 
a; = 0. Notation "b.t." stands for boundary terms which 
are irrelevant for the computation of the entropy. There- 



fore, the coefficients a^'"' and A^^' have the same geo- 
metrical form in this limit. 



{D) 
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The similarity of the heat kernel coefficients on smooth 
manifolds and manifolds with conical singularities at 
small angle deficits is a quite general property. If the 
components of the Riemann tensor are treated at the 
conical singularities as distributions |2l| the coincidence 
holds also for higher heat kernel coefficients and for dif- 
ferent field species j3^ . 

It should be noted that (|57|l does not hold for theo- 
ries with fields which have non-minimal couplings to the 
curvature. A typical example is the scalar theory with 
the wave operator A — ^i?, where ^ is a constant and 
R is the scalar curvature. The coefficient a[^'' does not 
depend on ^ because, by the construction of the effective 
action r(/3) in Section II. B, the wave operator on M^f^ 
coincides with the Laplacian A. However, the coupling 
changes the heat kernel coefficients on curved manifolds. 
In particular, the multiplier ^ in A\^^ in H53|) is replaced 
byi-e 

It is still an open question whether the entanglement 
entropy in a physical theory depends on the non- minimal 
couplings. The answer to it may be related to the pro- 
cedure of the measurements of the energy (a discussion 
of the topic can be found in In what follows we 

consider theories without non-minimal couplings. 



B. The gravitational coupling and analogue models 
of gravity 

It is natural to suggest that a fundamental gravity the- 
ory and a simple model with a cutoff considered in the 
previous Section share two features: i) in the both cases 
the low-energy gravity sector is a pure effective theory 
and equations for the metric tensor are determined en- 
tirely by the polarization properties of the physical vac- 
uum; ii) the theories do not have a problem of the ultra- 
violet divergences. 

The string theory which is considered as one of the can- 
didates to the role of the fundamental theory has these 
features. The second property ensures that the induced 
Newton coupling Gat is finite and can be expressed in 
terms of microscopical parameters. If the non-minimal 
couplings are absent Gn can be derived from the ground 
state entanglement in a flat spacetime by using H57|l . 

To get more intuition about quantum gravity it makes 
sense to study other models where the effective Newton 
coupling can be calculated with the help of (|57() . One 
of the possibilities is to consider condensed matter sys- 
tems which are higher-dimensional generalizations of spin 
chains discussed in Section II. A. The requirements to the 
models are: 1) the models are lattice theories; 2) they 
have critical points which correspond to second order 
phase transitions; 3) near the critical points the models 
are described by relativistic QFT's with massive fields. 

The cutoff in this type of models is the lattice spac- 
ing. The entanglement entropy is finite and it can be 
computed numerically. The third requirement guaran- 



tees that in the near critical regime one can use results 
of Section II. D to get the leading asymptotic of the en- 
tropy. One can use for this purpose an effective QFT (not 
the underlying lattice theory itself). Then the density of 
the entropy per unit area yields the induced gravitational 
couphng (relation ifFTI) '). 

One of examples of such systems are higher- 
dimensional Ising models which are known to be equiv- 
alent in the critical regime to scalar field theories with 
self- interactions 0| . The models are also interesting for 
other reasons. The behavior of the most part of known 
second-order phase transitions is equivalent at the critical 
point to a three-dimensional (3D) Ising model js^. There 
are also indications that Ising models can be represented 
as theories of random (hyper) surfaces. In particular, as 
was conjectured in 33] near the point of the second order 
phase transition the 3D Ising model might be equivalent 
to a non-critical fermionic string theory. 

The two-dimensional Ising model is exactly solvable. 
The reduced density matrix for the ground state entan- 
glement can be diagonalized. This property significantly 
simplifies computations of the entropy. In higher dimen- 
sions exactly solvable Ising models are not known. Thus, 
one has to rely on numerical methods. 

It should be noted that the role of condensed matter 
systems as gravity analogs is known for many years but 
in a different context: to study quantum effects in the 
external gravitational fields ^34] . One of such effects is 
the Hawking radiation. Under certain conditions sound 
waves in a liquid Helium or in Bose-Einstein condensates 
behave as scalar excitations propagating in an effective 
curved background with a metric similar to that near a 
black hole horizon (see, e.g., j35j and references therein). 

Let us emphasize that because we are interested only 
in the behavior of the Newton coupling we do not need 
to introduce any metric (or its analog) in the condensed 
matter system. To determine the coupling it is sufficient 
to study the response of the effective action in flat space 
to introduction of conical singularity. 



IV. APPLICATIONS 
A. RG-flow of the gravitational coupling 

The knowledge of entanglement entropy as a function 
of parameters of the theory can be used to find renor- 
malization group (RG) evolution of the induced Newton 
constant. This information is important for understand- 
ing the behaviour of gravitational interactions at different 
scales. 

Consider as an illustration a one-dimensional spin 
chain. It has a compact momentum space with the ra- 
dius p = 2TT / g determined by the lattice spacing g. The 
Wilson RG transformation implies integration over high- 
energy modes with momenta in the interval t~^p < p < p 
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(with t > 1) followed by rescaling, p to p' — tp. As a re- 
sult one gets a theory with larger masses. For the Ising 
model with A 7^ 1 this is equivalent to increasing 
the difference |A — 1|. The RG-transformation drives the 
theory away from ultraviolet fixed point A = 1. 



S 




1 ^ 

FIG. 4: This figure is based on results of [3^ . It shows the de- 
pendence of the entanglement entropy in the Ising spin chain 
@ at a fixed A'^ as a function of magnetic field strength A. 
The critical point A = f is ultaviolet fixed point. The arrows 
show directions of the RG flow from the ultraviolet to infrared 
regions. 



The RG-evolution of the entanglement entropy for the 
Ising model is known precisely |36j. The dependence of 
the entanglement entropy in the Ising spin chain at 
a fixed as a function of magnetic field strength A is 
shown on Fig. 01 This behaviour is in accord with a 
general property: in unitary theories the entanglement 
entropy should not increase along the RG flow because 
RG transformations eliminate contribution of the high 
energy modes. 

The fact that the entropy is not increasing does not 
imply the same property for the induced coupling G^^ 
defined by H57|l . The coupling behaves as the density of 
the entropy, G~^{t) = t^~^f{t), where D is the number 
of space-time dimensions and f{t) is a function which 
has the same RG-evolution as the entropy. In three- 
dimensional model discussed in Section II. D the induced 
gravitational coupling defined by (|57() increases under 
RG-evolution. On the other hand, in four-dimensional 
model the gravitational coupling decreases which means 
that gravitational interactions get weaker in the infrared 
region. 



B. The scaling hypothesis 

The scaling hypothesis in classical critical phenomena 
asserts that the physics is determined by large scale fluc- 
tuations which do not depend on the underlying micro- 
scopical details. For instance, in a magnetic at temper- 
ature T in an external field h which undergoes a second 
order phase transition the physical properties are deter- 
mined by the large domains of aligned spins. The mi- 
croscopic atomic scale does not enter thermodynamical 



relations near the critical point T = Tc- Non- analytic de- 
pendence of the physical quantities on jT — Td is entirely 
determined by a "singular part" of free energy density 
/sing. This part has a universal scaling 

where ^ ~ \T — Tc\^^ is the correlation length and dh is 
the scaling dimension of h. Functions f± are universal in 
a sense they coincide for different systems from the same 
class. 

Let us denote by iS,„„ the surface density of the " renor- 
malized" part of the entanglement entropy in the ground 
state, see Section II. E. This is the part of the entropy 
which remains finite in the limit g = 0. In general, S,^„ 
is a function of external fields h, and other parameters, 
g, driving the phase transition of the system at g — g^. 
By using analogy with the classical scaling one can con- 
jecture that near a critical point S,^„ has the following 
scaling [s^ 

SU9,h)^r^''-''>s±{hr'') , (59) 

where £, = \g — gc\~'^ is the correlation length, is the 
scaling dimension, and s± are universal functions. It can 
be also suggested that the divergent part of the entropy 
is analytic in g. 

One can check the validity of these statements in 
scalar field models discussed in Section II. D. In the limit 
when the edge size L of the cube is large one finds 
from (123), (El that s± = -f/(24) in D = 3 and 
s± = -i/(247r)ln(m//z) in D 4. 

In analogue models of gravity the universality conjec- 
ture H59|l means that the part of the induced Newton 
coupling (|57|l which does not depend on the cutoff is a 
universal function. This part of the coupling can be de- 
termined entirely in terms of an effective low-energy QFT 
and it does not depend on the details of the underlying 
microscopic theory. 

C. The problem of black hole entropy 

Studying entanglement in gravity analogs might be 
helpful for understanding the microscopical origin of the 
Bekenstein-Hawking entropy S^^ of black holes. As was 
discussed in Section LB, if the gravity is entirely induced 
by some underlying degrees of freedom the entropy of a 
black hole can be related to the entanglement between 
observable states and states hidden inside the horizon. 
Understanding the relation between the two entropies in 
the framework of a local relativistic quantum field the- 
ory is plagued by the problem of the ultraviolet diver- 
gences. The definition of the induced Newton coupling 
either requires introduction of the ultraviolet cutoff or 
working with a special class of ultraviolet finite theories 
with non-minimal couplings [T^ . |l3l |. The presence of 
non-minimal couplings makes statistical interpretation of 
S^" a difficuh task. 
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The equality of S to the fundamental entanglement 
entropy S is implied in suggestion |0J. One of the argu- 
ments why it should be true is that both in classical and 
in quantum theory the Bekenstein-Hawking entropy can 
be derived by the conical singularity method from the 
gravitational action (see e.g. 0, UM)^ ''^ the same way 
which we employed for calculation of S. 

If equality Q holds one has a tool to learn more about 
the degrees of freedom responsible for the origin of S^^ 
by using gravity analogs. There are different questions 
which can be addressed in different models. For instance, 
can one relate these degrees of freedom in the 3D Ising 
model with non-critical stings, do the near- horizon sym- 
metries play any role in the entropy counting as was sug- 
gested in )37l |. and other questions. 



V. SUMMARY 

This paper has two purposes. First, we wanted to 
show how the effective action approach and the meth- 
ods of spectral geometry can be applied for derivation of 
the entanglement entropy in condensed matter systems 
in the regimes when these systems can be described by 
relativistic QFT's. In particular, such regimes appear in 
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